Induced chiral Dirac fermions in graphene by a periodically modulated magnetic field 
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The eflect of a modulated magnetic field on the electronic structure of neutral graphene is exam- 
ined in this paper. It is found that application of a small staggered modulated magnetic field does 
not destroy the Dirac-cone structure of graphene and so preserves its 4-fold zero-energy degener- 
acy. The original Dirac points (DPs) are just shifted to other positions in k space. By varying the 
staggered field gradually, new DPs with exactly the same electron-hole crossing energy as that of 
the original DPs, are generated, and both the new and original DPs are moving continuously. Once 
two DPs are shifted to the same position, they annihilate each other and vanish. The process of 
generation and evolution of these DPs with the staggered field is found to have a very interesting 
patten, which is examined carefully. Generally, there exists a corresponding branch of anisotropic 
massless fermions for each pair of DPs, resulting in that each Landau level (LL) is still 4-fold de- 
generate except the zeroth LL which has a robust 47it-fold degeneracy with nt the number of pairs 
of DPs. As a result, the Hall conductivity a^y shows a step of size Antc^ /h across zero energy. 



PACS numbers: 73.43.Cd, 73.22.Pr, 73.61. Wp 

Low energy physics of neutral graphene is character- 
ized by the two inequivalent Dirac cones which is re- 
lated by the time-reversal symmetry and described by 
the relativistic massless Dirac equationjii^i^iii Nearly all 
important properties of neutral graphene is governed by 
the chiral massless fermions around the two cones. For 
example, the zero-energy anomaly due to the linear en- 
ergy dispersion and the particle-hole symmetry of the 
' Dirac cones give rise to the anomalous quantum Hall 
. effect (QHE)^iii^i^ii^ or the so-called half-integer QHE, 
I where the Hall conductivity is quantized to be half- 
integer multiples of Ae'^/h. When the Dirac-cone topol- 
ogy is destroyed or replaced by other structures, the sys- 
tem will undergo quantum phase transitions. In bilayer 
graphene, each Dirac cone is replaced by two touching 
I parabolic bands which leads to the 8-fold degen- 
eracy of zero-energy level;^ giving rise to the quantized 
Hall conductivity in bilayer graphene taken on values of 
integer multiples of Ae'^ /hMi^ 

Modulation of electronic structure in graphene has 
already been experimentally realized, where periodic 
electronioi^ii^ or magneti o^'^d^ potentials can be applied 
to graphene by making use of substrat o^^d'^d^dQ or con- 
] trolled adatom deposition;^ or by fabrication of periodic 
1 patterned gate electrodes. This kind of graphene super- 
lattice potential can change the Dirac-cone structure of 
graphene dramatically^!^ which may lead to some new 
phenomena, as well as potential application of graphene 
materials. 

In this paper, we present a study on the electronic 
structure of monolayer neutral graphene and its un- 
usual integer QHE under the influence of a periodically 
modulated orbital magnetic field, which is schematically 
shown in Fig.l. This kind of one-dimensional modula- 
tion of magnetic field can be achieved in experiments 
by applying an array of ferromagnetic stripes with al- 



ternative magnetization on the top of a graphene layer, 
or by making use of cold atoms in a honeycomb opti- 
cal lattice^i^ii^ or "artificial graphene" realized in a 
nanopatterned two-dimensional electron gasi^ Our anal- 
ysis shows that generally the Dirac-cone structure can 
not be smeared out by this time-reversal invariant mag- 
netic field. Similar to the cases of periodic electronic 
potential )^^'^^ new DPs will be generated with varying 
the amplitude of the field. What's remarkable and dif- 
ferent is that the newly generated DPs together with the 
original DPs will move and evolve in k space with the 
field. This leads to a series of quantum phase transi- 
tions with each phase characterized by its unusual integer 
QHE, which is expected to be observed by Hall measure- 
ments. 

We start with the tight-binding model on a honeycomb 
lattice in the presence of a perpendicular, periodically 
modulated orbital magnetic field. The Hamiltonian is 
given by, 

F=-i ^ e--^clc, +H.c., (1) 

where c| (q) is an electron creation (annihilation) oper- 
ator on site i, and < ij > denotes nearest- neighbor pairs 
of sites. Here the spin index is suppressed since we do 
not consider the Zeeman splitting. The magnetic flux 
per hexagon (the summation of a^- along the six bonds 
around a hexagon) is given by ^ = (f> ± S, where (j) 
measures the uniform magnetic flux whereas S the stag- 
gered modulated flux, both of which are in units of (po/iir 
with (po the flux quantum. Hereafter energy is measured 
in unit of the nearest-neighbor hopping integral t. 

To begin with, let us consider the effect of the two 
simplest types of staggered magnetic fields (SMFs), in 
order to extract the main physics behind graphene under 
the influence of a modulated orbital magnetic field. The 
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FIG. 1: (Color online) Illustration of the rectangular sample 
of graphene under periodically modulated magnetic fields, (a) 
and (b) represent two simplest SMFs, where each white and 
yellow (grey) hexagon has a fiux 5 and —S, respectively. The 
numbers 1, • • • ,4 represent the inequivalent atoms in a unit 
cell. Each arrow represents a phase shift sufi'ered by electrons 
when hopping along the direction, which is S/4 in case (a), 
and S/2 in case (b). (c) represents a long-period staggered 
fiux applied to graphene with lattice period 3aoLx, where ao 
is the lattice constant, (d) is a corresponding experimental 
layout in which there is an array of ferromagnetic stripes with 
alternative magnetization on the top of a graphene layer. 



configurations of the two types are schematically shown 
in Fig.[TJa) and (b), respectively, where proper gauge has 
been chosen for each case. 

We first show the evolution of DPs from an analytical 
calculation for SMF-I shown in Fig. [TJa). The tight- 
binding Hamiltonian in fc-spacc can be written as 



n = 
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where 7^ = — 2ie~*^ cos{^ky + |) and rjk = — te*' 
The Hamiltonian 7i determines the energy spectrum of 
electrons in graphene under SMF-I. The system has a pe- 
riodicity of 27r as a function of S due to gauge invariancc, 
so wc restrict 5 to range from to 2tt. The solution to 
the DPs can be easily obtained: the original DPs located 



0, cos(V3fcj^) 



clIj krf- 

and the newly generated DPs located at = i'l'/S, 



cos|, for < (5 < 47r/3, 



cos{^/3ky) 



cosf , for 27r/3 < (5 < 27r. 



An overall picture of the evolution of DPs in magnetic 
Brillouin zone(MBZ) under SMF-I is shown in Fig. H 
When (5 = 0, the original pair of DPs (red filled circles) 
are located at (0, ±27r/3\/3) [Fig. [DJa)]. As 6 increases, 
the two DPs move against each other along the ky di- 
rection, and eventually they reach the center of MBZ at 
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FIG. 2: (Color online) Schematic evolution of DPs in MBZ 
with increasing staggered fiux S under SMF-I. "•" and "•" 
represent the original DPs in pristine graphene and the in- 
duced (additional) DPs respectively, while the arrows rep- 
resent their moving directions. The black circles denote 
isotropic Dirac cones whereas the green (grey) ellipses denote 
anisotropic Dirac cones. The coordinates of the four corners 
of MBZ are (±7r/3, ±7^/^/3). 



S = 47r/3 [Fig.djf)] and thereafter disappear [Fig. HJg)]. 
On the other hand, right a,t S = 2tt/3, one additional 
pair of DPs arc induced simultaneously at the four cor- 
ners of MBZ (see the black filled circles in the figure). 
Note that only two of the four induced DPs are inequiva- 
lent and so only one pair contributes to the system, since 
the four induced DPs are all located at the boundary 
of MBZ. With increasing 6 the two newly induced DPs 
move along the lines kx = ±7r/3 towards the edge cen- 
ter respectively [Fig. [Hc)-(h)]. When 5 = 2tt, the two 
DPs arrive at (±7r/3, ±27r/3-\/3). In this case, electrons 
hopping along the arrows shown in Fig. [TJa) will suffer 
an additional phase 7r/2, which is a pure gauge. Thus 
the system corresponding to this value of 6 [Fig. [2jh)] is 
actually physical equivalent to that of Fig. [l^a), because 
they differ only by a gauge transformation. 

It is shown in Fig. [2] that DPs not only annihilate in 
pairs but also emerge in pairs.— This is interpreted by 
the fact that the two DPs in each pair are connected to 
each other by the time reversal symmetry which is still 
preserved by the SMF. Now we lay out our numerical 
results to support these findings. We show for different 
d the energy dispersion near zero energy along the two 
lines kx = [Fig. [3i;a)-(e)] and k^ = ±7r/3 [Fig. [Sljf)- 
(j)] where DPs reside. Note that 6 — 27r/3 and 47r/3 
are two critical values at which the pair of induced DPs 
emerge and the original DPs completely superpose each 
other, respectively. Apart from the two critical values, 
as 6 increases from to 27r, the number of DPs changes 
from one pair to two pairs and then back to one pair. 
This interesting evolution of DPs will dramatically affect 
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FIG. 3: (Color online) (a)-(e) Electron energy near the DPs 
under SMF-I versus ky with kx = for various values of 
staggered flux S. The dashed line represents the boundary of 
MBZ. (f)-(j)The same as (a)-(e) but with = n/3. (k) Hall 
conductivity a^y under SMF-I, with (f) = 27r/768 for several 
values of S. Inset: The two renormalized factors j3 {f3') as 
functions of staggered flux S. 



the degeneracy of the LLs which can be reflected by the 
Hall conductivity. 

The Hall conductivity can be calculated directly 
through the standard Kubo formula^I by numerical di- 
agonalization of the Hamiltonian (1). In Fig. E^k), the 
resulting Hall conductivity a^y near zero energy is plot- 
ted as a function of the Fermi energy Ep. According to 
the Hall plateaus steps in axy, the system can be classi- 
fied into three types. For < S < 2tt/3 or 47r/3 < S <2Tr, 
with spin degeneracy taken into account, axy has a step of 
size Ae^/h, which is the same as that of pristine graphene, 
shown in Fig. [Stk) for S = OAn. For 2tt/3 < 6 < 47r/3, 
axy has a step of size 8e^//i across zero energy (neutral 
filling) whereas a step of size Ae^/h in the other energy 
range, which can be seen in Fig. [3Uk) for 5 = Q.Stt. Re- 
markably, right at (5 = TT, numerical results of Hall con- 
ductivity show that all the steps have the same size of 
8e'^/h. We interpret these phenomena as follows. 

The isotropic Dirac cones become anisotropic under 
the influence of a modulated magnetic field (see Fig. [5]). 
The chiral fermions around an anisotropic Dirac cone 
can be physically described by the anisotropic pseudospin 
Hamiltonian, 
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where p± = apx ± ibpy, vp = Mao/2h is the Fermi ve- 



FIG. 4: (Color online) Schematic evolution of DPs in MBZ 
with increasing staggered flux S under SMF-II. All symbols 
used here have the same meanings as that in Fig. [2] The 
coordinates of the four corners of MBZ are (±7r/3,±7r/2V3). 



locity, and the two dimensionless coefficients a and b 
measure the degree of anisotropy of the cone. In the 
presence of a uniform magnetic field B this anisotropy 
gives rise to a renormalized LLs En = iphvpy^^Jnl/lB 
with /3 = \/ab a dimensionless renormalization fac- 
tor, and Ib — \/4'o/ 47ri? the magnetic length. With 
spin degeneracy taken into account, it is found that 
for < S < 47r/3, the 4- fold degenerate LL spectrum 
near the original Dirac cones have a P value given by 
P - m = {73T+^(mf^)H^ (note when 
5 = 0, /5 = 1), while for 27r/3 < S < 2t:, the 4-fold 
degenerate LL spectrum near the induced Dirac cones 
have another different value given by (3' = /3(27r — 6). 
For a general 6 between 27r/3 and 47r/3, the LLs for the 
two branches are not degenerate except the zeroth LL, 
which is exactly 8-fold degenerate. The zeroth LL is in- 
dependent of the external uniform magnetic field so its 
8-fold degeneracy cannot be removed, leading to a %e? /h 
Hall conductivity step at the zeroth LL and a 4e^//i step 
at other LLs. However, when 5 = tt, the two factors 
are equal to each other, i.e., /3 = /?', all LLs for these 
cones (which are isotropic now) overlap and so are ex- 
actly 8-fold degenerate. Therefore at 5 ~ n., the Hall 
conductivity can be expressed as axy = 8{N + l/2)e^//i, 
with N LL index. We remark that actually, within the 
range of S where the two pairs of Dirac cones coexist, 
there should exist a series of critical values of 6 given by 
0^/13''^ = p/q, where p and q are two coprime integers. 
At these critical values, besides the zeroth LL, the mgth 
LL (with m = 0, 1, 2, . . . ) in the original pair of cones are 
exactly degenerate with the mpth LL in the induced pair 
of cones, giving rise to a 8e^//i Hall conductivity step at 
these energies. Another critical value of (5 is at (5 = 47r/3, 
where the zero point a,t kx = ky = Q is not a DP, but 
rather a semi-DP. Around the semi-DP, energy disper- 
sion is found to be linear along kx^ but parabolic along 
ky. This peculiar feature can be compared with the elec- 
tronic structure in V02-Ti02 nanoheterostructuresi^ 

Now we turn to explore the second type of the SMF 
shown in Fig. [TJb). Fig. 2] shows the schematic picture 
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FIG. 5: (Color online) (a)-(d) Electron energy near the DPs 
under SMF-II, as a function of ky with = Q for various 
values of staggered flux &. The dashed lines represent the 
boundary of MBZ. (e)-(h)The same as (a)-(d) but as a func- 
tion of kx, with ky = 7r/6\/3, instead, (i) Hall conductivity 
(Jiy under SMF-II with </> = 27r/768 for several values of 5. 



of the evolution of DPs in the corresponding MBZ. Like 
that of SMF-I, at the beginning of varying (5, the two DPs 
of pristine graphene are located at (0, ±7r/3-\/3), and then 
they move towards the origin along ky direction. When 
6 — 7r/2. each DP of the pair changes into three DPs 
at k==(0,±7r/6%/3) [Fig. [li;b)], but they completely su- 
perpose each other and can not be distinguished there. 
After that the original DPs go on moving along ky direc- 
tion and eventually arrive at the origin at 5 = 27r/3 and 
then vanish, whereas the other two pairs of induced DPs 
move along the ky = ztTr/OVS direction until 6 ^ n they 
reach the points (±7r/6, ±7r/6\/3) [Fig. [DJf)], and then 
backtrack. The evolution of DPs from 5 = tt to 27r is just 
the reverse of the above process. 

Compared with the SMF-I, in a period of 6, there are 
four critical values {S = 7r/2, 27r/3, 47r/3, and 37r/2), at 
which the DPs emerge or vanish. All the induced DPs 
under SMF-I move parallel to the ky axis while that un- 
der SMF-II move parallel to the axis. This should 
be associated with the configuration of the SMF, where 
the induced DPs incline to move towards the periodic 
direction of the SMF. 

The electronic energy spectrum near the zero energy 
are shown in Fig. [n)Ja)-(h). What is significant is that as 
i5 increasing from to tt, the number of DPs changes 
from one pair to three pairs and then to two pairs. 
So, the DPs indeed emerge and annihilate in pairs. In 
Fig. O^i) , the Hall conductivity cr^y is plotted as a func- 



TABLE I: The evolution properties of DPs in MBZ at < 
S < 2n/Lx for various Lx's {L^ is from 2 to teens). 
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"rio represents the number of pairs of DPs in MBZ with kx = 0. 
''rii represents the number of pairs of DPs in MBZ with kx = 
■izir/SLx. 

'^nt (= no + rii) represents the total number of pairs of DPs in 
MBZ. 

''<5c is a critical value where the number of DPs at the line fc^ = 
changes to zero. For Lx = 2,3, • ■ ■ , 10, &c is approximately equal 
to O.SSStt, 0.5957r, 0.4657r, 0.3757r,0.3157r, 0.2757r, 0.2457r, 0.2177r, 
0.1967r, respectively. 

> 2. For L is an odd number, rii = 0, 2, 4, ■ ■ . , L — 1, L; for L 
is an even number, nj = 0, 2. 4, . . . , L. 



tion of the Fermi energy Ep. For < 5 < 7r/2 or 
37r/2 < i5 < 27r, the Hall conductivity can be expressed 
as (Txy — 'i{N + l/2)e^//i, which is the same as that of 
pristine graphene. This means that the Dirac-cone topol- 
ogy is preserved within this range without new induced 
DPs. For tt/2 < S < 27r/3 or 47r/3 < S < 37r/2, a^y has 
a 12e^/h step across the zeroth LL and Ae^/h or 8e^/ft. 
step at the other LLs, implying the zeroth LL is 12-fold 
degenerate. Interestingly, for 27r/3 < 5 < 47r/3, all LLs 
are 8-fold degenerate and the Hall conductivity can be 
expressed as Uxy = 8(A^ -I- l/2)e^//i. This expression 
with quantized values of half-integer multiples of jh 
is robust and is interpreted by the fact that the original 
pair of DPs has disappeared and the induced two pairs 
of DPs are located symmetrically in MBZ giving rise to 
exact 8-fold degeneracy of their corresponding LLs. 

Thus far we have discussed two simplest SMFs where 
the magnetic flux alternates along the armchair chains or 
zigzag chains, respectively. Now we generalize our theory 
to the cases of SMF with long spatial period, shown in 
Fig. [TJc). Taking SMF-I for example, we make the mag- 
netic flux alternate every L^. zigzag chains. Numerical 
analysis shows that, as 8 increases from to 2tt /L^r^ 
the number of pairs of original DPs Uq (at k^ — 0) de- 
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FIG. 6: (Color online) The case with Lx = 6. Electron energy 
near the DPs under a long-period SMF shown in Fig.[T](c) , as 
a function of ky with (a)-(d) kx = and (e)-(f) kx = tt/18 for 
various values of staggered flux 5. The figures in (b),(f) and 
(c),(g) have been shifted to the left by 0.1447r and 0.0937r, 
respectively. The total number of pairs of DPs is 1,3,5,6, 
respectively for the four 5 values. 



creases from one to zero, while that of the induced DPs rii 
(exactly located at the boundary of MBZ = ±Tr/3Lx) 
increases from zero to gradually obeying a sequence 
of 0, 2, 4, • • • ,Lx. Accordingly, the total number of pairs 
of DPs is rit = 1, 3, 5, • • • , La; + 1, L^- For the detail, see 
Table I. In Fig. [6l we take = 6 for example. When 
< 6 < 2it/Lx, the number of pairs of the induced DPs 
shows a sequence of Ui = 0,2,4,6 [Fig. [6l^c)-(h)], while 
that of the original pair is always no = 1 or [Fig. [H^a)- 



(d)]. Once again, the DPs emerge and annihilate in pairs. 

After application of such a long-period SMF, each pair 
of Dirac cones has generally different anisotropy, i.e., has 
different renormalization factor /? in their corresponding 
LL spectrum, so all LLs except zeroth LL are still 4-fold 
degenerate. However, the zeroth LL is exact 4nt-fold de- 
generate, leading to a step of size Ante^ /h in a^y at the 
zeroth LL)^ as 5 increasing from to 2tt/ L^- Therefore, 
under a modulated orbital magnetic field, the property 
of neutral graphenc is actually governed by the number 
of pairs of DPs. This is a significant signature and can be 
detected by Hall measurements. So long as is no more 
than the magnetic length of the systems, the physics con- 
tained in them is similar. But for graphene system with 
the magnetic length much less than L^, the DPs will be- 
come more and more dense and will finally be merged 
into the zeroth LL of graphenci^ This deserve further 
study and will be discussed elsewhere. 

In summary, we have investigated the electronic struc- 
ture in neutral graphene under periodically modulated 
magnetic fields. It is found the modulated magnetic field 
can induce additional DPs in graphene and the evolution 
of these DPs can be manipulated by the magnitude and 
period of the field. These induced DPs add additional 
degeneracy to the LLs, especially a 4nt-fold degeneracy 
at the zeroth LL, leading to an unusual integer QHE near 
neutral filling. These phenomena are expected to be ob- 
served by Hall measurements. 
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